AO-A09Q  337 
UNCLASSIFIED 


SHEFFIELD  UN IV  <ENGLANO)  DEPT  OF  CHEMICAL  ENGINEER IN — ETC  F/8  *1/5 
A  THREE-DIMENSIONAL  MODEL  OF  SPRAY  EVAPORATION  IN  GAS  TURBINE  C— CTC(U) 
APR  80  W  H  AYERS  AFOSR-80-017* 

HIC-343  AFOSR-TR-80-092B  NL 


AFOSR.TR-  80  -  0^  2  8 


urn 


*  £ 

*  Vv  .»  r--. 

o.  v,  3 1:  \ 


& 


■/ •  r?1 


\ 

/  X"  ! 


/lF0tf-W'  &7<f 


'V' 


/,.  A  Three-Dimensional  Model  of 

'is  j 

_y  Spray  Evaporation  in  Gas 
Turbi neCombu s t  o  r s • 


W#  H, 'Ayers ^ 

Aprrtpr^®80^ 


WJWfcr  -Ow  >  <  r 

<3l:|tr  4  *-,<  •  ,  f  j  .  ,  . 


80  9  22  ll3f 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  or  This  PACE  (Whan  Data  Entered) 

REPORT  DOCUMENTATION  PAGE 


[  1.  REPORT  NUMBER 


AFOSR.TR.  80-092 


J  DArP  head  instructions 

*  rMUu _ _ _ rr.F CRT  COMPLETING  FORM 

/]2.  GOV'T  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 

/  1  AC'tiC'9#337 


|«.  TITLE  (and  Subtitle) 


A  THREE  DIMENSIONAL  MODEL  OF  SPRAY  EVAPORATION  IN 
GAS  TURBINE  COMBUSTORS 

7  Au7hOR?7)  ~  "  — -» ~  • 

W  H  AYERS 


5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

Interim 


6.  PERFORMING  ORG.  REPORT  NUMBER 

_ HIC  343  _ 

8.  CONTRACT  OR  GRANT  NUMBERS 


AF0SR-80-0174  7 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  10.  PROGRAM  ELEMENT,  PROJrCT.  TASK 

AREA  ft  WORK  UNIT  NUMBERS 

UNIVERSITY  OF  SHEFFIELD  '' 

MAPIN  STREET  2308/A2 

SHEFFIELD.  Si  3JD _ 61102F _ 

11.  CONTROLLING  OFFICE  NAME  ANO  AOORESS  12  REPORT  DATE 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH/NM  April  1980 

BUILDING  410  13  number  of  pages 

BOLLING  AFB.  DC  20352 _ 33 _ 

U  “monitoring  AGENCY  NAME  ft  ADDRESS^//  different  from  Controlling  Office)  15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

1 Sa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

17  distribution  statement  (of  this  Report ]  " 


2308/A2 
61I02F _ 

12.  REPORT  DATE 

April  1980 

13.  number  OF  PAGES 


Approved  for  public  release;  distribution  unlimited. 


*7.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 


IB.  supplementary  NOTES 


[  19.  KEY  WORDS  (Confirm*  on  reverse  side  if  necessary  and  identify  by  block  number ) 


SPRAY  GAS  TURBINE  COMBUSTION 

COMBUSTION  TOO  PHASE  FLOW 

DROPLETS 

THJIEE  DIMENSIONAL  COMPUTATIONS 

r MATHEMATICAL  MODELING _ 

BSTRACT  /'ConUnue  on  reverse  aide  it  necessary  and  identify  bv  block  numM' 

The  result  of  this  study  is  a  computer  program  which  can  calculate  the  trajectory 
of  a  droplet  or  particle  in  three  dimensions  in  a  three  dimensional  flow  field. 
Droplets  are  injected  into  a  computer  predicted  flow  field  in  a  gas  turbine 
combustor  can  with  initial  conditions  corresponding  to  the  point  of  sheet  break 
up  of  fuel  leaving  the  atomizer  being  modeled.  The  fundamental  equations  of 
motion  are  solved  numerically  for  each  droplet  size  in  a  statistical  distribu¬ 
tion  representing  the  whole  spray.  The  angular  position  oi  the  point  ol 

injection  is  also  varied.  The  rate  of  fuel  evaporation  in  each  cell  of  the 


DO  ,:r73  1473 


_ UNCLASSIFIED  / 

SECURITY  CLASSIFICATION  of  THIS  PAGE  -Wh*n  nm«  l  ntered 


% 


Summary 


The  result  of  this  study  is  a  computer  programme  which  can 
c-lculate  the  trajectory  ofadroplet  or  particle  in  three  dimensions 
in  a  three  dimensional  flow  field* 

Droplets  are  injected  into  a  computer  predicted  flow  field 
in  a  fas  turbine  combustor  can  with  initial  conditions  correspond} ng 
to  the  point  of  sheet  break-up  of  fuel  leaving  the  atomiser  being 
modelled.  The  fundamental  equations  of  motion  are  solved  numerically 
for  each  droplet  size  in  a  statistical  distribution  representing  the 
whole  spray.  The  angular  position  of  the  point  of  injection  is  also 
varied. 

The  rate  of  fuel  evaporation  in  each  cell  of  the  finite  difference 
grid  is  then  calculated  so  that  this  can  be  used  for  improved  flow  field 
calculations.  Droplet  trajectories  are  also  presented  graphically 
to  give  information  on  the  three-dimensional  nature  of  the  spray. 
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Notation 


cfc 

°o 

cr 

Jr 
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Nu 

Pr 

Re 

t 

T 
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n,  v, 
u,  v, 

y> 
•  * 
X,  y, 


Evaporation  constant 
Drag  coefficient 

Specific  heat  capacity  at  constant  ^pressure 
Droplet  or  particle  diameter 
Rosin-Rammler  mean  diameter 
Latent  heat  of  evaporation 
Mass 

Mass  flowrate 

Rosin-Rammler  spread  parameter 

Nusselt  number 

Prandtl  number 

Reynolds  number 

Time 

Temperature 

Resultant  velocity  vector 

w  velocity  components  in  axial,  radial  and  tangential  directions* 

w  time  derivatives  of  u,  v  and  w 

z  distance  in  axial,  radial  and  angular  directions* 
z  time  derivatives  of  x,  y  and  z* 


P 


Density 

Viscosity 

Thermal  conduct ivity 


Subscripts: 

O  Initial 

to  Flov;  field 
^  Gas 

(  Liquid 

p  Particle  or  droplet 
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Introduction 


The  availability  of  large  general  purpose  digital  computers 
has  enabled  numerical  techniques  to  be  applied  to  complex  physical 
phenomena#  In  the  field  of  Combustion  Research,  mathematical 
modelling  can  be  used  to  predict  the  behaviour  of  all  kinds  of 
combustion  processes  and  can  give  information  on  the  chemical  and 
physical  processes  taking  place  (Ref#  l)#  As  a  design  tool, 
mathematical  modelling  can  assist  in  providing  a  less  empirical 
approach  to  the  problems  of  combustor  design,  a  field  of  increasing 
importance  as  a  result  of  the  present  interest  in  improved  fuel 
efficiency  and  reduction  of  pollutant  emissions. 

The  object  of  this  study  is  to  produce  a  general  algorithm 
to  calculate  the  trajectory  of  a  particle  in  a  gas  flow  in  three 
dimensions#  The  equations  of  motion  and  trajectory  are  solved 
numerically,  as  an  analytical  solution  is  not  practicable  in  three 
dimensions  with  the  coordinate  system  used* 

The  algorithm  is  then  applied  to  the  specific  case  of 
evaporating  fuel  droplets  in  a  gas  turbine  combustor,  in  order  to 
model  mathematically  the  spray  evaporation  taking  place  and  to 
provide  information  on  the  trajectories  and  behaviour  of  individual 
droplets  and  size  groups  with  varying  initial  conditions# 

In  order  to  model  the  evaporating  droplets,  equations  representing 
droplet  heating  up  and  evaporation  must  be  included  in  the  trajectory 
algorithm# 

The  computer  predicted  gas  flow  field  for  a  Ivcoming  combustor, 
used  in  the  study,  was  produced  at  Sheffield  (ref.  ?),  and  defines 
the  geometry  of  the  finite  difference  grid  into  which  droplets  are 
injected# 

Work  in  this  field  has  already  produced  droplet  trajectories 
in  two  dimensions  (Ref#  5)* 


M 
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Having  obtained  trajectories  of  individual  droplets,  the 
model  is  extended  to  an  entire  spray  by  calculating  the  trajectories 
of  a  range  of  droplet  sizes  and  initial  conditions.  The  trajectories 
are  presented  graphically  to  yield  information  on  the  behaviour  of 

i 

individual  droplets  within  the  spray. 

The  rate  of  fuel  evaporation  in  each  grid  cell  can  also  be 
calculated  for  the  fuel  spray#  This  information  can  th eh  be  used 
in  the  original  flow  field  calculations,  replacing  the  time  consuming 
method  currently  used,  i.e#  representing  each  droplet  size  range  ns  a 
chemical  species  and  solving  the  resulting  elliptic  equations. 


-  /  - 

Theory* 

Because  of  the  cylindrical  nature  of  combustion  chambers 
and  the  gas  flows  involved,  a  cylindrical  polar  coordinate  system 
is  used  to  describe  the  location  of  the  particle  or  droplet,  where 
x  is  the  axial  distance  from  the  datum,  y  is  the  radial  position, 
and  z  (radians)  is  the  angular  position*  Velocities  are  measured 
in  the  axial  (u), radial  (v)  and  tangential  (w)  directions. 

The  equations  of  motion  of  a  particle  ,  neglecting  all  forces 
except  drag,  in  component  form  are: 


dXp  -  —  F  (  ^  p  ~  ^ 


where  ^p  Xf,vr  are  the  absolute  particle  velocities  in  the  axial, 
radial  and  tangential  directions  respectively,  and  f 

are  the  corresponding  gas  velocities. 

For  the  general  case,  F  is  given  by: 


j 


where  Dp  is  the  particle  diameter,  He  is  the  relative  Reynold's 


number,  defined  ass 


Re  =  hh. 


where  -Up 

and  are  the  resultant  velocities  of  the  particle  and 

gas  stream  respectively. 

The  drag  coefficient,  Cp  ,  is  given  by  the  following  equations 
(  Ref.  3)  : 

Cd  ~  27  Re  O^Re^>$0(6) 


Co  =  0-271  fa. 


O'Zii 


80<Rc^lO^  (7) 


Co  -  2 


10^<  Re. 


If  Stoke*  s  regime  prevails  and  there  is  uniform  gas  flow 
(  i.e.  are  constant  ),  then  the  axial  component  can 

be  solved  analytically. 


unity 


f  Co  Rt\ 

Under  Stoke*  s  regime,  the  term  I - I  is  approximated  to 

\  24  / 


Assuming  D p  remains  constant,  the  term 


( |  B  \ 


becomes 


constant,  whereupon  equation  (l)  is  integrable. 


In  a  Cartesian  coordinate  system,  the  equations  of  motion 
for  all  three  components  would  have  the  form  of  equation  (l),  thus 
permitting  an  analytical  solution  for  obtaining  the  droplet  or 
particle  trajectories#  Hov/ever,  in  the  rectangular  polar  coordinat 
system  used,  the  interdependence  in  the  expressions  for  the  redial 
and  tangential  components  prevents  an  analytical  solution  being  found 
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For  the  rectangular  polar  coordinate  ‘.ystem,  the  equations 
of  trajectory  are: 


(9) 


<jr  =  \  (10) 

= 

Equations  (l)  to  (ll)  describe  the  trajectory  of  a  particle 
in  a  gas  stream.  For  an  evaporating  droplet,  Df  is  not  constant, 
and  the  rate  of  change  of  droplet  diameter  vrith  time  is  required. 

For  forced  convection  this  ir  given  by: 


Pp  _  _  b 

dt  ~  JPp 

v:here  is  the  evaporation  constant.  This  is  dependant  on 

the  properties  of  the  fuel  car-  vrell  as  the  surrounding  gases, 
and  is  given  by  ’dice  .and  Agoston  as  (  Ref.  4)  : 


(1+  0-23  Re'")  (12) 


Pt  q*} 


(13) 


If 


(l  +  O  25 


he  a 


can  be  integrated  to  give: 


slimed  constant, 


then  equation  (l2) 


Pf!  =  Dp<?  -  C/,(l+  O'ZSRe'-’Xt-t.)  (i/!) 


I 
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ropre renting  the  droplet  r1i' meter  Pp  at  time  t,  ’here  Pp0  i ^ 
the  droplet  diameter  -it  time  t0  . 


.  To  allow  for  droplets  "hooting  up"  a  ft  or  entry  into  the  hot 
p:as  stream  it  has  been  assumed  tha,t  droplets  do  not  be£dn  to  evaporate 
until  they  reach  their  boiling  point.  The  rate  of  chnnye  of 
temperature  v:ith  time  is  £iven  hy: 


He  c  6  /Va-V  .  ( Xo  -  TP) 

it  Pi  v?  c'< 


(15) 


vjhere  Tp  is  the  droplet  temperature.  Assuming  (  _  tL t  _  1  1 

V  Pc  Pp  / 


is  constant,  on  integration  this  gives; 


T,  =  Xo-  iXo-Trc) 


<l6) 


vrhere^Ipp  is  the  drop  temperature  at  time  t0  . 


Equation  (14)  is  suppressed  until  the  droplet  temperature, 
reaches  the  boiling  point  of  the  liquid#  The  husselt  number,  IJu 
is  given  by: 


Nii  r  Z  ^  0*G  Re  2  Fr  *  (17) 


v:here  the  Prandtl  number,  Pr  is  given  by: 

Ph  -  f!l  (. 

X? 

The  above  solutions  are,  of  course,  valid  only  when 
(t-t0 )  is  taken  to  be  sufficiently  small. 


\ 
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Application* 


Equations  (l)  to  ( 1 8 )  above  describe  the  motion  of  an 
evaporating  droplet  in  three  dimensions  ih  a  gas  stream. 

Numerical  integration  can  now  be  used  to  calculate  the  droplet 
trajectory  in  the  given  gas  flow  field. 

For  the  purposes  of  this  invest igation  a  computer 
predicted  flow  field  in  a  gas  turbine  combustor  was  used*  This 
has  been  calculated  for  a  Lycoming  combustor  using  a  numerical 
prediction  algorithm  to  model  the  physical  and  chemical  processes 
taking  place  (  Ref*  2  ). 

The  geometry  of  the  Lycoming  combustor  can,  which  was  designed 
as  a  research  combustor,  is  shown  in  Fig,  1.  Air  enters  via  a 
swirler  which  produces  a  finite  swirl  velocity*  The  fuel  spray  is 
introduced  at  a  point  near  the  centre  line  in  front  of  a  haffle 
and  is  of  hollow  cone  type.  Primary,  secondary  and  dilution 

air  streams  are  introduced  by  three  sets  of  injection  holes,  each 

of 

consisting 'six  equally  spaced  circular  orifices  around  the  periphery 
of  the  combustor  can.  This  is  responsible  for  the  three-dimensional 
nature  of  the  problem. 

For  the  computer  model,  the  combustor  can  is  represented 
by  a  single  60°  sector  containing  one  set  of  air  injection  holes 
to  represent  the  six  identical  segments.  This  is  then  subdivided 
into  a  further  7  sectors*  In  the  axial  direction  the  length  of 
the  combustor  is  divided  into  27  sections,  and  radially  into  18 
sections,  giving  a  total  of  3402  elements  in  the  finite  difference 
grid.  The  grid  is  shown  in  Fig.  2* 

The  elements  are  labelled  (  I,J,K  )  in  the  x,y  and  z 
directions.  Note  that  in  the  angular  direction  elements  with 
K  between  2  and  6  inclusive  occupy  10°  sectors  whereas  the  K=1 
and  K=7  elements  occupy  a  sector  of  5°  •  In  fact  the  two  5° 
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sectors  are  adjacent  and  can  be  considered  as  one  10°  sector, 
but  due  to  the  positioning  of  the  grid  and  in  order  to  simplify 
calculations  they  are  considered  as  two  separate  half  sectors. 

For  the  purposes  of  this  investigation  the  flov;  parameters 
used  (  which  are  the  three  components  of  velocity,  temperature 
and  density)  are  defined  for  each  grid  element  and  are  assumed 
to  be  valid  throughout  this  control  volume.  The  complex  flow 
throughout  the  combustor  can  thus  be  represented  by  the  finite 
difference  grid. 

The  flow  data  is  obtained  from  the  output  of  the  combustor 
model  described  in  Ref.  2.  The  same  grid  geometry  is  used  in 
the  latter  but  the  gas  velocities  are  defined  for  a  staggered 
grid.  This  is  demonstrated  by  Fig.  3(a)  which  shows  that  the 
u  ,  v  and  w  velocity  control  volumes  do  not  coincide  with  those 
for  temperature,  pressure  etc.  In  order  to  considerably  simplify 
the  representation  of  velocity  data  the  u,  v,  w  velocity  control 
volumes  have  been  redefined  to  coincide  with  the  nominal  control 
volumes  by  averaging  adjacent  velocity  values,  an  example  being 
shown  in  Fig.  3(b). 

Since  the  first  set  of  air  injection  holes  occurs  at  the 
primary  air  inlet  it  might  be  expected  that  the  gas  flov:  field  would 
be  essentially  two  dimensional  upstream  from  this  point.  It  is  on 
this  assumption  that  the  two  dimensional  model  of  Ref.  5  is 
justified.  To  investigate  the  validity  of  this  assumption  a 
statistical  analysis  of  u-,  v-  and  w-velocity,  temperature  and 
density  data  was  made.  For  each  value  of  I  and  each  of  three 
J  values,  the  set  of  data  for  all  K  values  was  averaged  and  a  sample 
standard  deviation  calculated. 


Sample  standard  deviation,  s,  is  given  by: 


(19) 


n  - 1 
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where  x,  ,  xz . x*  are  the  values  of  the  data,  x  is  the  mean 

of  the  data  and  n  is  the  number  of  data  points.  ;’ince  we  are 

interested  in  the  amount  of  angular  variation  as  a  proportion  of 

the  mean  of  the  values  rather  than  the  absolute  variation,  the 

sample  standard  deviation  is  expressed  as  a  percentage  of  the  mean. 

\ 

Hence: 


Samp 1 e  St  anda rd 
Deviation  (Percent) 


(20) 


The  results  of  this  analysis  are  shown  in  Figs.  4(a)  to  4(e) > 
and  are  discussed  below. 


A  vector  plot  of  the  velocity  distribution  is  also  given  in 
Pigs.  5(a)  and  5(^)*  It  can  be  clearly  seen  that  there  is  a 
recirculation  zone  in  front  of  the  baffle,  and  that  the  flow  in 
the  region  of  the  air  injection  holes  is  dominated  by  the  inv:ard 
flow  of  air.  The  flow  cross-section  of  Pig.  5(b)  at  the  point 
of  primary  air  inlet  clearly  shows  how  the  three-dimensional 
nature  of  the  flow  field  is  brought  about  by  the  disturbance 
introduced  by  the  air  injection. 


Once  the  u-,  v-,  w-velocity,  temperature  and  density  data 
is  stored  as  a  27  x  18  x  7  array,  the  fundamental  equations  of 
motion  (  equations  (l),  (2)  and  (3)  )  can  be  solved  numerically 
by  assuming  that  the  gas  velocity,  temperature  and  density  remain 
constant  throughout  each  cell. 

These  three  simultaneous  differential  equations  are  solved 
by  using  the  fourth  order  Runge-Kutta  method.  The  time  step 
length  is  calculated  to  give  approximately  the  same  number  of 
time  steps  while  the  droplet  is  in  each  cell.  After  each  time 
step  the  equations  of  trajectory  (  equations  (9)?  (10)  and  (ll)  ) 
are  used  in  finite  difference  form  to  calculate  the  new  values  of 
x,  y  and  z  as  follows: 


Xft  -  X K-,  -f  (  U.K  +-  tin-,) 


(21) 


I 
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In  addition,  for  an  evaporating  droplet,  either  equation 
(l4)  or  (l6)  is  used,  depending  on  whether  the  droplet  has  reached 
its  "boiling  point,  to  calculate  droplet  diameter  or  temperature. 

Any  auxiliary  equations  vdth  terms  containing  u,  v,  w,  x,  y  or  7. 
or  Df  are  also  calculated  after  each  time  step. 

A  check  is  also  made,  after  each  time  step,  to  establish 
whether  the  particle  or  droplet  has  passed  through  the  boundary 
of  the  volume  cell.  If  this  proves  to  be  the  case  the  new  cell 
is  found  and  new  values  for  the  constants  for  the  cell  are  obtained* 

Also,  Ppo  of  equation  (14)  and  Tfo  of  equation  (l6)  and  the  corresponding 
t0  are  reset.  A  new  time  step  length  is  calculated  to  maintain  an 
approximately  constant  number  of  time  steps  in  each  cell*  A  check 
is  also  made  for  the  following  specific  cases: 

1)  Droplet  passes  into  an  adjacent  60°  sector# 

In  this  case  the  z  value  of  the  droplet  is  adjusted  to 
bring  it  to  the  corresponding  position  in  the  original 
sector,  i.e.  60°  is  subtracted  from  its  z  value* 

The  actual  angular  position  is  retained  for  output* 

2)  Droplet  passes  through  centre-line. 

Although  there  would  seem  to  be  a  very  small  chance  of 
this  occurring  a  provision  is  incorporated  for  the  droplet 
to  be  reflected  from  the  centre-line  if  the  radial  distance,  y, 
should  become  less  than  zero* 

3)  Droplet  hit 3  a  wall* 

Since  the  original  velocity  data  is  not  valid  for  cells 
for  which  J  is  greater  than  17*  a  droplet  is  said  to 
have  hit  the  combustor  wall  if  it  enters  a  cell  for  which 
J  is  equal  to  l8* 
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4)  Droplet  passes  through  the  exit  of  the  combustion  chamber. 

If  a  droplet  or  particle  survives  to  pass  through  a  cell 

for  which  I  =  ?7  it  is  reported  as  having  left  th**  combustion 
chamber • 

\ 

5)  Droplet  evaporates. 

If  the  diameter  of  an  evaporating  droplet  falls  to  1  cap 
than  (lO)^  metres  (  approximately  3  microns  )  it  is 
reported  as  having  evaporated.  In  practice,  droplets 
vanish  almost  immediately  as  Pp  becomes  very  srmll,  and 
droplets  smaller  than  this  behave  erratically  unless  the 
time  step  length  is  reduced  greatly.  This  test  is  made 
with  each  time  step. 

6)  Droplet  remains  in  combustion  chamber  after  5^0  time  steps. 

The  iteration  is  halted  if  more  than  500  time  steps  are 
required  in  order  to  prevent  an  excessive  use  of  computer 
time . 


Any  of  conditions  3)  to  6)  above  will  terminate  the  iteration 
loop,  otherwise  the  next  time  step  is  calculated. 


When  a  droplet  is  evaporating  the  amount  of  fuel  evaporated 
per  unit  time  in  each  cell  is  calculated  as  it  leaves  the  cell, 
given  bys 


(?4) 


■where  Diiv  =  droplet  diameter  on  entry  into  cell 
=  droplet  diameter  on  exit  from  cell 
58  droplet  diameter  on  start  of  trajectory 
=  total  mass  flowrate  of  fuel  represented  by  the 
droplet,  e.g.  the  mass  flowrate  of  fuel  in  the  rise  range  being 
represented  by  a  single  droplet. 


-  13  - 

If  the  droplet  hit?  a  wall  as  in  3)  above,  the  rerr  ining 

fuel  contained  in  the  droplet  is  "evaporated11  into  the  cell 

adjacent  to  the  wall,  representing  droplets  adhering  to  the  wall 

of  the  combustion  chamber#  It  would  be  possible  to  extend  the 

model  to  include  droplets  rebounding  off  the  wall,  shattering  etc# 

* 

To  extend  the  model  to  predict  the  behaviour  of  a  whole  spray, 
a  statistical  size  distribution  is  used  to  represent  the  spray  ar  a 
number  of  size  ranges,  each  represented  in  turn  by  a  single  droplet 
diameter,  the  mean  of  the  size  range. 

It  is  assumed  that  after  sheet  break-up,  the  fuel  spray 
leaving  the  atomiser  obeys  a  Rosin-Rammler  size  distribution. 

In  this  model,  the  mass  fraction  of  fuel  having  a  droplet  diameter 
greater  than  D  is  given  as: 

=  ~'HiJ }  (J5> 

where  D  is  the  Rosin-Rammler  mean  diameter  and  n  is  the  spread 

parameter#  For  this  study  a  value  50  microns  was  chosen  for  5 
and  the  spread  parameter  is  taken  as  2.  The  mass  fraction  of 
fuel  in  a  given  size  range  is  found  by  subtraction: 

"w  -|  e^Ht)7  - exf HtJ  I  <?6) 

vrhere  the  size  range  lies  betvreen  D,  and  Dj  • 

The  size  ranges  are  chosen  to  be  of  equal  length  and  cover 
droplet  diameters  from  zero  to  a  diameter  belov;  which  995^  of  “the 
mass  of  fuel  falls#  A  total  of  20  size  ranges  was  chosen  to 
adequately  model  the  complete  spray# 
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oirt^le  droplets  with  diameters  equal  to  the  ru  an  of  each  size 
ranf-e  are  injected  into  the  fpis  flow  field  with  initial  x  and  y 
coordinates  corresponding  to  the  point  of  sheet  break-up  of  tVu  spr  ay. 

In  addition,  for  each  size  ranp*e ,  droplets  '  re  injected  from  si  r 
different  angular  locations  differing  by  10  .  Consequently  the 
trajectories  of  a  total  of  1?0  fuel  droplets  are  calculated. 

The  atomiser  characteristics  are  specified  by  spray  cone  an  pi  o  ,  v*  1  oc 
distance  from  nossle  to  sheet  break-up,  and  the  Rosin-Rami  er  me-  n 
and  spread  parameter#  From  this  and  pener  1  physical  d  la  the 
entire  spray  cone  is  constructed. 

The  amount  of  fuel  evaporated  in  each  cell  ir  summed  for  Ihe 
whole  spray.  The  mass  fl  or ’rate  of  fuel  represent  rd  by  e*  rh  ■"o^let 
trajectory  is: 


(r) 


where  in  defined  in  eruption  Since  e-,ch  ri  re  r  nff 

represented  by  six  droplets  of  different  initial  nr. pul  nr  position, 


factor  of  1  /6 


The  total  mass  flowrate  of  fuel  is 


i  s 
ihe 


occurs 
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Results  and  Pi  scusri  on. 

Th^  ~erultr  of  th°  rt '  t  i  rt  i  c*'l  'r.^yrir  of  4h^  ~ng,'!~r 
v-rintion  t  ho  flow  ^.rmot^rr  ir  presented  in  Fi  ■ — .  /!  (  n  ) 

to  d(r).  Percentage  -tand-^d  devi~tionh  r> ro  given  for  J  valuer 
of  3*  0  and  16  nr>  indicated  in  the  key  (  Fig.  /?(f)  ).  It 
can  be  seen  that  angular  variation  of  all  five  par*. net  err  in 
generally  loss  than  20jS  upstream  of  the  primary  ^ir  injection, 
at  which  point  it  begins  to  increase  very  rapidly.  It  ^ould 
therefore  appear  that  in  the  region  downstream  of  the  grid 
location  at  v:hich  I  ir>  equal  to  10  the  flov:  field  can  he 
considered  to  he  three-dimensional  and  to  he  two-dimensional 
upstream  from  this  point. 

The  results  of  calculations  for  fuel  sprays  of  evaporating 
droplets  are  shO’Tt  graphically  in  Figures  5  to  9  as  follows: - 

inclirled  cone  angle  spray  velocity 


( degrees) 

(m/s) 

PiC*  5 

45 

20 

PiC*  6 

45 

10 

Fig.  7 

45 

50 

PiC*  8 

30 

20 

PiC*  9 

8o 

20 

Distance  from  atomiser  nozzle  to  sheet  hreak-up  is  taken  as 
5mm.  The  number  of  size  ranges  is  20  except  for  Fig.  8  where 
it  was  reduced  to  12  to  reduce  the  amount  of  computer  time  required. 
Other  relevant  data  can  he  found  in  the  programme  listing  of  Appendix  1* 

The  results  are  shown  as  an  orthographic  projection  of  the 
droplet  trajectories  consisting  of  a  side  elevation  of  the  combustor 
and  an  end  elevation  looking  along  the  axis  of  the  combustion  chamber. 

The  air  injection  locations  are  indicated  by  arrows.  Droplet 
burn-out  locations  are  encircled.  Droplets  whose  loci  do  not  end 
in  a  circle  have  hit  the  combustor  wall. 


1 
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The  combustor  dimensions  are  shown  in  Fig.  1*  Part  of  the  section 
of  the  combustor  can  downstream  of  the  primary  air  injection  holes 
is  not  sho;.7i  in  Figs*  5  to  9  • 

Figure  5  shows  the  predicted  rprav  obtained  for  an  included 
nominal  cone  angle  of  45°  anl  initial  sprnv  velocity  of  TO  m/*', 
these  characteristics  being  typical  of  the  pres  sure- jet  type  of 
atomiser*  It  can  be  seen  that  a  large  number  of  droplets 
penetrate  into  the  three-dimensional  region.  Those  droplets 
passing  within  the  vicinity  of  the  primary  air  stream  are  deflected 
and  some  pass  right  across  the  combustion  chamber  and  impinge  on 
the  opposite  wall.  Similar  behaviour  is  observed  in  Fig,  6  where 
the  spray  velocity  is  only  10  m/s  ,  as  a  result  of  which  the 
droplets  have  evaporated  after  travelling  a  shorter  distance  along 
the  chamber.  Figure  7  shows  a  higher  velocity  spray  (  5 0  m/s  ) 
which  has  resulted  in  a  spray  which  remains  in  liquid  form  further 
along  the  combustion  chamber ,  some  droplets  passing  beyond  the 

secondary  air  injection  holes. 

Figure  8  shows  the  behaviour  of  a  spray  with  a  narrower 
cone  angle.  The  spray  has  been  broken  up  to  a  greater  extent 
than  for  the  45°  cone  angle  because  the  majority  of  the  droplets 
have  entered  the  recirculation  zone  where  there  are  greater 
velocity  gradients*  It  can  be  seen  that  some  droplets  show  an 
apparently  erratic  behaviour  as  they  are  repeatedly  deflected 
by  several  primary  and  secondary  air  streams. 

Figure  9  shows  a  spray  with  a  large  cone  angle  (  80°), 

In  this  case  no  droplets  have  penetrated  into  the  three-dimensional 
7, one  and  a  large  proportion  of  the  droplets  have  collided  with  the 
combustor  wall  upstream  of  the  primary  air  injection  holes. 

Figure  10  is  a  histogram  shrvinr  Ihe  mass  flow  rate  of 
fuel  (area  on  graph)  against  distance  along  the  x-xin  of  the 
combustor*  This  has  been  calculated  for  the  3^°  cone  angle  case 
of  Fig,  8,  It  shows  that  a  significant  proportion  of  Ihc  fu^l 
evaporates  in  the  three-dimensional  region  of  the  conbu  *  or  can. 
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In  view  of  the  fact  that  come  droplets  show  a  very  abrupt 
deviation  from  their  original  path  in  the  vicinity  of  the  air 
injection  streams,  a  check  was  considered  necessary  to  ensure  that 
this  was  not  due  to  inaccuracies  in  the  solution  procedure  arising 
from  too  large  a  time  step  length*  The  time  step  length  was 
therefore  reduced  by  a  factor  of  four,  and  several  sets  of 
trajectory  calculations  were  performed  to  verify  the  original 
results*  There  was  no  significant  change  in  the  results* 


I ; 


i 


Conclusions 


It  has  been  demonstrated  that,  according  to  the  computer- 
predicted  flow  data,  the  conditions  in  t,he  region  of  the  primary 
air  injection  holes  and  downstream  from  this  point  are  three-dimensional 
in  nature.  It  has  also  been  shown  that  it  is  possible  to  produce  a 
three-dimensional  mathematical  model  of  individual  droplet  trajectories 
within  a  fuel  spray  and  thus  to  model  the  spray  as  a  whole.  In  the 
absence  of  any  experimental  data  it  has  not  been  possible  to  perform 
a  comparison  betv:een  the  predicted  and  actual  observed  spray  behaviour. 
Experimental  data  to  verify  the  predictions  is  still  required* 

It  can  he  seen,  however,  that  the  techniques  developed 
could  be  of  use  in  combustor  design.  It  is  clearly  undesirable 
to  have  droplets  impinging  on  comhustor  wells  causing  • wet1  spots, 
or  to  have  un-evaporat ed  fuel  surviving  as  far  as  the  lower  temperature 
regions  of  the  combustor  where  it  may  lead  to  incomplete  combustion 
and  consequently  poor  fuel  efficiency  and  the  presence  of  pollutants 
in  the  exhaust  gases*  The  effect  of  cone  angle  and  spray  velocity 
on  these  phenomena  has  been  shown*  Further  study  along  these  lines 
might  enable  a  more  suitable  atomiser  design  to  be  found  or  lead  to 
improvements  in  combustion  chamber  design  with  particular  regard  to 
the  positioning  of  the  air  injection  ports. 

It  must  be  noted  that  although  the  droplet/particle  trajectory 
algorithm  has  been  applied  to  the  specific  case  of  evaporating 
droplets  in  gas  turbine  combustors,  it  remains  a  general  procedure 
and  is  capable  of  being  used  to  mode]  the  motion  of  any  particle 
or  droplet  in  a  three-dimensional  flow  field.  Examples  of  applications 
include  diesel  engines  and  cyclone  separators* 
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1)  Geometry  of  Lycoming  combustor  can. 

2)  Finite  difference  grid  for  the  Lycoming  combustor. 

3)  (a)  Location  of  u,  v  and  w-velocity  control  volumes. 

(b)  Relocation  of  u-velocity  control  volumes. 

4)  (a)  Angular  variation  of  u-velocity. 

(b)  Angular  variation  of  v-velocity. 

(c)  Angular  variation  of  w-velocity. 

(d)  Angular  variation  of  temperature. 

(e)  Angular  variation  of  density. 

(f)  Key. 

( g )  Velocity  Field  :  Position  K  =  4* 

(h)  Velocity  Field  :  Fosition  I  =  13* 

5)  Droplet  trajectories  for  spray:  cone  angle  =  45°*  velocity  =  20  m/s. 

6)  Droplet  trajectories  for  spray:  cone  angle  =  45°>  velocity  =  10  rn/s. 

7)  Droplet  trajectories  for  spray:  cone  angle  =  45°>  velocity  =  50  m/s. 

8)  Droplet  trajectories  for  spray:  cone  angle  =  30°*  velocity  =  20  m/s. 

9)  Droplet  trajectories  for  spray:  cone  angle  =  80°,  velocity  =  20  m/s. 

10)  Histogram  of  mass  flow-rate  of  fuel  vapour  vs.  distance  along  x-axis. 
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